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^.^ . Abstract 

CO , Here I discuss the major pitfalls and the most severe mistakes of the above mentioned paper. 

The thorough analysis shows that despite all the claims the present work has nothing to do with 
the thermalization process in relativistic heavy ion collisions. In contrast to the authors' beliefs 
I show that their main result is not derived, but is a combination of mathematical mistakes and 
fj . hand waving arguments. 
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Thermalization phenomenon observed in heavy ion cohisions at intermediate and high energies 
is one of the pillars of modern nuclear physics. It means a transformation of two colliding pieces of 
cold nuclear matter, which produce the highly excited non-equilibrated strongly interacting matter, 
into a media being in local thermal equilibrium. There exist a large number of theoretical works 
devoted to its explanation on a basis of various approaches and models. The latter are varying 
from the gluon transport model [1] and hadronic cascade model [2, 3] to rather sophisticated hydro- 
kinetic approaches [4]. 

In the present work [5] the authors claimed no less than they found "an explanation for the 
thermalization phenomenon, starting from such basic concepts as the conservation laws of energy 
and momentum and leaving the details which are characteristic of every specific process, aside for 
further researches". However, their model, the maximal isotropization m,odel (MIM), has nothing to 
do with the collisions of heavy ions and, as shown below, does not shed any light on a thermalization 
process at all. 

The authors, indeed, account for the energy and momentum conservation in all two particle 
collisions and in this way they redistribute the total energy and momentum of the system among 
the ideal Boltzmann particles. However, the authors have neglected any inelastic reaction in the 
system and production of any particle in it. Therefore, from the beginning the possible range of 
applications, if we consider the hadrons, is restricted from above by the lab. energy per particle 
which must be below the pion mass m.,^ ~ 140 MeV. However, the authors of [5] continue to use 
the relativistic expression for particle energy and even apply their model to pions with the mean 
energy up to 800 MeV. After I realized this, it became clear to me that the thorough discussion of 
such 'tiny issues' as spherical symmetry of the momentum distribution after each elastic collision 
assumed by authors, or an absence of any cross-sections in all their equations, or the energy and 
time dependence of the momentum integration volume Vp, introduced by them, e.t.c, would require 
too much time, and hence I decided to analyze just a few most important issues. 

The authors operate with the simple equations which are well known to the community [6, 7], 
but fail to analyze or even simply to discuss the validity of their main approximation made in 
Eq. (15). However, the approximation suggested by the authors drastically changes the properties 
of the system from the normal ones to the anomalous ones. It is well known from the standard 
course of statistical mechanics (see, for instance, paragraph 67 of [8]) that, if the microcanonical 
density of states increases or decreases as a certain power of particle energy, then the system has 
the normal properties since the corresponding canonical system can have positive temperatures 
only. Unfortunately, the microcanonical density of states obtained after the approximation (15) for 
fixed values of a„ and 6„ is exponentially decreasing function of particle energy and hence it leads 
to the anomalous system properties since it allows the existence of negative temperatures [8]. This 
mistake has the fatal consequences for the present paper. 

However, the wrong consequences of the approximation (15) are not seen immediately because 
the authors make another mistake by loosing their original object of analysis. Thus, originally 
the authors are dealing with the microcanonical distribution (11) and then suddenly they make 
the Laplace transform to the canonical density of states in Eq. (12). The rest of the paper is 
devoted to the analysis of the canonical distribution (12) with the temperature of the external 
thermostat given by the parameter of the Laplace transform T = P~^ which afterwards appears in 
all relevant expressions. Thus, instead of elucidating the temperature of the microcanonical system 
which would be necessary for studying the thermalization phenomenon, the authors insert this very 
temperature T = j3~^ by hands through the canonical partition. Therefore, the present ■work 
does not give any novel piece of correct information on the thermalization process in 
heavy ion collisions. 

Perhaps, the authors of MIM do not know that, when the Laplace transform (12) is made, 
it is assumed that there exists an ensemble of microcanonical systems of all allowed energies and 
each of these systems is brought into a contact with the external thermostat whose temperature 



is T = I3~^ (for a detailed discussion see [9, 10]). It seems that the authors of this work do not 
know that the microcanonical ensemble allows one to define the microcanonical temperature Tmce 
via the derivative of the entropy Smce with respect to the total energy Etot in a standard way as 
'^mce = ( dE'"^ I [9> 10]- If they applied such a formula to the microcanonical density of states, 
it would not be necessary to write the wrong and meaningless discussion given in Sect. 3.4 (see 
below) . 

Moreover, the authors do not understand that to elucidate the temperature of the system of 
colliding particles one has to analyze the microcanonical density of states, since in the canonical 
ensemble the outer thermostat forces system to adopt its temperature. The Laplace transform 
to the canonical ensemble is an auxiliary mathematical trick to simplify the calculation of the 
original microcanonical partition with the constraints (the conservation laws in this work). After 
the constraints are evaluated in the canonical ensemble, one has to perform the inverse Laplace 
transform to the original ensemble and analyze it. Some general results on the both Laplace 
transforms in the semiinfinite or finite interval can be found in Refs. [11] and [12], respectively. I 
am sure that these references would be very useful to the authors of [5], since they do not know 
that in the right hand side of the inverse Laplace transform there must be an additional factor g^j^, 
which is missing in their Eqs. (52), (53) and in the discussion presented between Eqs. (51) and 
(53) (see Sect. 3.4.). This mistake was extremely surprising to me since one of the authors of Ref. 
[5] was present at many of seminars of mine where the results of works [11] and [12] just devoted 
to various applications of the Laplace transform technique were discussed in great details! 

Furthermore, the authors of [5] claim that they derived the nonequilibrium momentum dis- 
tribution function. The same statement is repeated in their subsequent work [13]. As was shown 
above none of the distributions listed in Ref. [5], was actually derived. Moreover, by claiming 
this the authors clearly indicate an absence of any knowledge of great amount of works devoted 
to studies of different aspects of nonequilibration phenomena occurring in relativistic high energy 
nuclear collisions. The authors of Ref. [5] implicitly mean that the 'derived' nonequilibrium mo- 
mentum distribution function has some correction to the Boltzmann momentum distribution of 
relativistic particles of ideal gas. Note, however, that more realistic interaction between relativistic 
particles (reltivistic hard-core repulsion) leads to a far more complicated momentum distribution 
[14] than the Boltzmann one and this fact has nothing to do with the collision of two relativistic 
nuclei, but is an inherent property of the equation of state of hadronic matter. Thus, accounting 
for realistic interaction between hadrons one obtains the modification of the Boltzmann momentum 
distribution. It is also well known that the non-Boltzmann momentum distribution of the secondary 
hadrons formed in the high energy nuclear collisions can be a consequences of a complicated col- 
lective motion of particles [15] which have the Boltzmann distribution in their local rest frame, 
but are measured by detector in a lab. frame. In addition, the conditions of a kinetic freeze-out 
process of the secondaries also affect their momentum spectra [16, 17, 18, 19, 4] and this modifies 
the Boltzmann distribution of particles in the fluid elements freezing out at the time-like freeze-out 
hypersurfaces [17, 18]. 

From the discussion above it is evident, if the authors of [5] applied the standard thermodynamic 
relation to the microcanonical density of states obtained from the canonical one after the inverse 
Laplace transform over f3 with exp[(3Etot], then they would get the result independent of the Laplace 
transformation parameter P = T^^ . Hence I conclude that the main result of the present paper does 
not correspond to the collisions of heavy ions at all, but it represents the mixing of two (directed to 
each other) slow flows of small hard balls in a box. The box and balls are in a thermal contact with 
an external thermostat while the interaction with the box's walls is neglected. Thus, the present 
work does not corresponds to its title and to the claims written in its introduction. Moreover, it 
scientific value seems to be negative. 

As explained above, the MIM does not account for particle production and, hence, it cannot 
be used to relativistic nuclear collisions. Nevertheless, without any explaining remark the authors 



apply their analysis to the ideal (Boltzmann) gas of pions at ultrarelativistic energies per particle 
(up to 800 MeV). However, besides other mistakes, this implies that the heavy ions consist of free 
pions. 

All this is not surprising because, in contrast to the MIM, in high energy nuclear collisions a 
thermalization occurs not only due to elastic scattering of the particles, but also and mainly due to 
production of new (sorts of) particles [20, 3] which automatically increases the local particle density 
in a system which, in its turn, essentially enhances the reaction rate leading to a thermalization. 
Since such a production of particles is not and cannot be accounted within the framework of the 
present work, its output is unphysical. 

Also I have to stress a very sloppy style of presentation and its mentor tone which are in a 
tremendous dissonance with the trivial mistakes and pitfalls of this work. Thus, one hand the 
present authors pretend on a solving a fundamental problem, but on the other hand they are not 
careful with their definitions, approximations and even wording. After I read a few paragraphs 
and did not find any word of what kind of statistics is used, I was not already surprised by such 
pearls as "the confinement of the momentum space" or "the nonnormalized distribution" (see the 
sentence after Eq. (10)) and so on. 

In summary, the most original part of present work is a collection of primitive mistakes which 
have nothing to do with the thermalization problem in heavy ion collisions. 
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